Dynamics of a dislocation bypassing an impenetrable precipitate: the Hirsch 

mechanism revisited 



Takahiro Hatano 

Earthquake Research Institute, University of Tokyo, Tokyo 113-0032, Japan. 
(Dated: February 2, 2008) 

Dynamical process where an edge dislocation in fee copper bypasses an impenetrable precipitate 
is investigated by means of molecular dynamics simulation. A mechanism which is quite different 
from the Orowan mechanism is observed, where a dislocation leaves two prismatic loops near a 
precipitate: i.e. the Hirsch mechanism. It is found that the critical stress for the Hirsch mechanism 
is almost the same as the Orowan stress, while the spatial inhomogeneity of the shear stress is 
essential to the Hirsch mechanism. We also find that the repetition of the Hirsch mechanism does 
not increase the critical stress. 
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Plastic flow in crystalline materials is mainly domi- 
nated by the movement of dislocations. They are curvi- 
linear defects which can move at much lower stress lev- 
els than the theoretical strength of a perfect crystal Q. 
A dislocation interacts with other lattice defects such 
as voids or precipitates, which obstruct the dislocation 
motion to result in hardening. Such dislocation-obstacle 
interactions dominate plasticity of crystalline materials, 
and hence they are one of the major concerns in mate- 
rials science. Until very recently, they were considered 
exclusively by plausible continuum models @, , which 
deal with a single glide plane and neglect atomistic dis- 
creteness. The recent computer development enables di- 
rect molecular dynamics (MD) simulation on dislocation 
systems, which consist of multi- million atoms. MD sim- 
ulations have revealed dynamical properties and atom- 
istic details of dislocations. For example, edge disloca- 
tions absorb vacancies when they interact with voids 
or stacking fault tetrahedra 0. They are good illustra- 
tions of the usefulness of MD simulation in the field of 
dislocation physics. Along the line of these studies, the 
interaction with an impenetrable precipitate, which has 
not been investigated by MD simulation, is explored in 
this Letter. 

An interaction between a dislocation and a precipitate 
is characterized by their shear moduli. When the shear 
modulus of a precipitate is larger than the bulk's, the 
interaction is repulsive so that the precipitate can be im- 
penetrable 6]. In this case, the extent of hardening be- 
comes significant. This is the principle of so-called "parti- 
cle strengthening" , which is widely utilized in processing 
stronger materials 7]. There is a mechanism by which 
a dislocation bypasses impenetrable obstacles, where a 
dislocation largely bows out to leave a dislocation loop 
around the obstacle |l|. Note that the loop is on the 
original glide plane. This process and the resultant dis- 
location loop are referred to as the Orowan mechanism 
and an Orowan loop, respectively. However, because the 
interaction between an Orowan loop and a dislocation 
is strongly repulsive, theoretical calculations in which 
Orowan loops are assumed to be persistent predict un- 
reasonably strong work-hardening H, . To avoid this 



contradiction, some alternative bypass mechanisms have 
been presented. Hirsch 0, had pointed out the pos- 
sibility of prismatic loops formation in a bypass mecha- 
nism. This mechanism is named after Hirsch and referred 
to as the Hirsch mechanism. On the other hand, Brown 
and Stobbs presented another mechanism where a dislo- 
cation loop of a secondary Burgers vector nucleates on 
the precipitate surface [l^- In both cases, the resultant 
structure is a row of prismatic loops which was exper- 
imentally observed by transmission electron microscopy 
(TEM) [HI, E3, However, since the dynamical pro- 
cess cannot be observed by TEM, the postulated mech- 
anisms have not been directly tested. In this Letter, we 
wish to clarify the dynamics of prismatic loops formation 
and to find out quantitative conditions which determine 
resultant dislocation loops: Orowan loops or prismatic 
loops. Along the line of thought, an MD simulation is 
performed. 

Let us describe the computational model. We consider 
fee copper utilizing a many-body interatomic potential 
of Finnis- Sinclair [l5| and adopt parameters which are 
determined by Ackland et al 16]. The lattice constant 
a = 3.615 A. The dimensions of the model system are 
x = 23 x [112] (20 nm), y = 79 x [110] (40 nm), and 
z = 27 x [111] (17 nm). This system consists of approx- 
imately 1.2 million atoms. (In addition to this system, 
we prepare some smaller systems which have different 
length in the x direction.) We introduce an edge disloca- 
tion parallel to the x axis, setting the Burgers vector to 
be parallel to the y axis. Note that we focus on an edge 
dislocation, since a screw dislocation can avoid a precipi- 
tate by switching its glide planes via cross- slip P3, which 
is not our interest here. An impenetrable precipitate is 
modeled by a set of immobile atoms, which are coherent 
to the bulk crystal. It mimics a precipitate of the infinite 
shear modulus, which cannot be sheared by dislocations. 
Periodic boundary conditions are employed in the x and 
the y directions, and the surfaces exist only in the ±z di- 
rections. In order to cause the shear strain, (111) surface 
is displaced at a constant velocity. We test two cases, 
which are referred to as conditions I and II, respectively. 
In condition I, the upper (+z) surface move towards the 
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FIG. 1: Successive snapshots of an edge dislocation bypass- 
ing an impenetrable precipitate. The radius of the precip- 
itate (blue atoms) is 1.5 nm. A dislocation (red atoms) is 
visualized by omitting atoms which have 12 nearest neigh- 
bors, (a) A dislocation bows out to form a screw dipole. (b) 
The screw dipole cross-slips, (c) The screw dipole undergoes 
double cross-slip to annihilate each other at the top of the 
precipitate, (d) The dislocation has bypassed the precipitate 
with a superjog. A prismatic loop is left behind. 
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FIG. 2: Geometry of the primary cross-slip from (111) to 
(111), (a) A [111] projection, (b) A [111] projection. 



—y direction while the lower surface is fixed. In condi- 
tion II, both surfaces move opposite to each other. In 
both conditions, we set the strain rate e = 7 x 10 6 [s -1 ]. 
Note that we adopt the condition I unless explicitly in- 
dicated. Temperature is set to be 300 K, where the ve- 
locities of copper atoms are given randomly according to 
the Maxwell-Boltzmann distribution. 

Let us turn to the simulation results. First, we focus 
on the dynamics and the geometry. The configuration is 
shown in Fig. ffl 17], where a dislocation does not leave 
an Orowan loop but two prismatic loops. At the first 
stage, a dislocation bows out to form a screw dipole of 
opposite signs. The main difference between the Orowan 
mechanism and the mechanism observed here lies in the 
behavior of this screw dipole. It cross-slips onto a (111) 
plane (see FIG.[2|), and then undergoes double cross-slip 
onto another (111) plane, on which the screw dipole even- 
tually annihilates each other. In FIG.Q(d), a prismatic 
loop is formed on the left side of the precipitate and a su- 
perjog is formed on the opposite side. Then the superjog 
is dragged by the dislocation for a while, and is eventu- 
ally separated to form the secondary prismatic loop as 
shown in Fig. [31(a). If the primary cross-slip takes place 
towards the extra atomistic half plane which constitutes 
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FIG. 3: (a) Configuration just after a single bypass process. 
A dislocation leaves two prismatic loops: an interstitial loop 
and a vacancy loop. The secondary loop is formed relatively 
far from the precipitate, since it has been dragged by the 
dislocation as a superjog. (b) A row of prismatic loops after 
the passage of three dislocations. Note that the loops on the 
right side of the precipitate do not form an apparent row. 



the edge dislocation, the prismatic loop created before 
the precipitate is interstitial and the secondary loop af- 
ter the precipitate is vacancy. (If the primary cross-slip 
goes opposite, the order is reversed). Note that the re- 
sultant prismatic loops have the same Burgers vector as 
that of the original edge dislocation. Indeed, this mecha- 
nism is identical to the Hirsch mechanism [ll| . We stress 
that this is the first case where the Hirsch mechanism is 
dynamically observed, while TEM observations have de- 
tected only the resultant structures. Note that the mech- 
anism of Brown and Stobbs has not been observed 
in the conditions investigated here. More importantly, 
it is found that the Hirsch mechanism always occurs for 
the condition I, while the Orowan mechanism is realized 
for the condition II. To clarify the underlying physics, we 
have to investigate the quantitative aspect. 

We investigate the critical resolved shear stress for the 
Hirsch mechanism, which we will call the " Hirsch stress" . 
The shear stress is defined as the area-averaged force 
acting on the surfaces. The Hirsch stress is defined as 
the maximum shear stress during a bypass process. It 
is realized just before the primary cross-slip, by which 
the Hirsch stress is determined. Behaviors of the Hirsch 
stress with respect to the precipitate radius and to the 
precipitate spacing are shown in Fig. 01 where we find 
the Hirsch stress is described by the following logarithmic 
law. 
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where A and B denote material-dependent parameters, 
r is the radius of a precipitate, and L is the precipitate 
spacing in the x direction. Namely, L = L x — 2r, where 
L x denotes the system size in the x dimension. Note 
that A = 6.5 N/m and B = 0.5 nm in the present sys- 
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FIG. 5: Spatial profiles of the average atomistic displacement 
in the y direction when the shear stress becomes maximum. 
The red solid line represents 5y(z) for the Hirsch mechanism, 
while the green dashed line denotes that for the Orowan mech- 
anism. The gradient of Sy(z) for the Hirsch mechanism is 
asymmetric with respect to z = 95 A (the glide plane). Sys- 
tem parameters r = 1.5 nm and L = 20.4 nm. 



FIG. 4: The critical resolved shear stresses for the Hirsch 
mechanism (denoted by the red x symbols). Note that the 
Orowan stress is almost the same as the Hirsch stress. The 
green dashed lines represent Eq. Q with A = 6.5 N/m and 
B = 0.5 nm. (a) Dependence of the critical stress on the 
precipitate radius, where L x is 20.4 nm. (b) Dependence of 
the critical stress on the spacing of the precipitates, where r 
is 1.5 nm. 



tern. This kind of logarithmic behavior of the critical 
stress is universal in the context of dislocation-obstacle 
interactions 0, ^| , the reason of which is explained as 
follows. From the balance of forces that are transverse 
to the dislocation line, the resolved shear stress (r) and 
the radius of curvature (p) of a bowing dislocation are 
interrelated as 
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Note that 7 denotes the line tension, which depends on 
the configuration of a bending dislocation. When an edge 
dislocation bends to form a screw dipole at a pinning 
point, the line tension is estimated as [2| 
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Since a fully bending dislocation forms a semi-arc be- 
tween two pinning points, the radius of curvature is 
p = L/2. Inserting this relation and Eq. ((SJ into Eq. 
immediately leads to the logarithmic law. 
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Although Eqs. and Q are of the same form, 

the quantitative discrepancy is considerable. Namely, 
Gb/2ir = 1.96 N/m in Eq. while the correspond- 

ing factor in Eq. Q is A = 6.5 N/m, which is three 



times larger than Gb/2ir. There are at least two reasons 
for this disagreement. The first is that dissociation of dis- 
locations is not taken into account in deriving Eq. It 
takes much more stress to bend two partial dislocations 
simultaneously pj|. The second is the high dislocation 
density. Since the glide of a single dislocation causes large 
strain relaxation in small systems, the critical depinning 
stress may depend on the initial position of a dislocation 
in MD simulations. 

Then we discuss the difference between the Hirsch and 
the Orowan mechanisms, both of which are realized in 
the present simulation depending on the boundary con- 
ditions. It should be remarked that the Orowan stress 
is almost the same as the Hirsch stress, and therefore 
we must refine the analysis. Recalling that the Hirsch 
mechanism is realized by displacing only the upper sur- 
face, inhomogeneity of the shear strain may play a key 
role. In order to see the inhomogeneity more explicitly, 
we calculate the average atomistic displacements in the 
y direction for each mechanism. We define 
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where the subscript i denotes the i-th atom. We take the 
sum regarding all the atoms that belong to the layer of 
z < Zi(t) < z-\- 1.89 A. Then n(z) is the number of atoms 
in the layer, and t is the time just before each bypass 
mechanism begins (i.e. when the shear stress becomes 
maximum.) FIG. shows Sy (z) for each mechanism. For 
the Hirsch mechanism, the gradient of Sy is different with 
respect to the glide plane. Because the cross-slip which 
precedes the Hirsch mechanism occurs towards the larger 
stress region, it is concluded that the strong spatial in- 
homogeneity of the shear stress is essential to the Hirsch 
mechanism. This spatial inhomogeneity is attributed to 
the existence of the precipitate. Because the precipitate 
consists of immobile atoms, it inhibits the deformation 
caused by the boundary displacement. Therefore, if the 
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upper (or lower) surface is displaced, the elastic strain 
tends to localize on that side [2l[ . When the both surface 
is equivalently displaced, the elastic strain becomes sym- 
metric to realize the Orowan mechanism. The above dis- 
cussion further leads us to the speculation that Orowan 
loops become unstable for the larger shear strain. They 
can be decomposed into two prismatic loops via double 
cross-slip of the screw component, which is essentially the 
same manner as the Hirsch mechanism. For example, in 
the system shown in FIG.0 the shear strain of 0.103 is 
enough to destabilize Orowan loops. 

While the discussions so far involve a single bypass 
process, for the rest of this Letter, iteration of the Hirsch 
mechanism is briefly investigated. Namely, we consider 
the dynamics of subsequent dislocations which succes- 
sively interact with a precipitate. It involves the rele- 
vance of the Hirsch mechanism to work-hardening. Al- 
though prismatic loops can move under an appropriate 
stress field (prismatic punching), they cannot move in 
the simple shear condition to be stagnant near the pre- 
cipitate. However, they do not cause hardening since 
subsequent dislocations can easily penetrate them. For 
example, when L x = 20 .4 nm and r = 1.5 nm, it takes 
less than 180 MPa for dislocations to penetrate prismatic 
loops, which is much smaller than the Hirsch stress itself, 
560 MPa. Plastic flow is maintained by repetition of the 
Hirsch mechanism, which yields a row of prismatic loops 
as is shown in Fig. H9(b). This is what has been observed 
in the experiments [Tj, [Tj, [l3( . By contrast, the loops on 
the right side of the precipitate do not form an apparent 
row. Because the prismatic loops move towards the +y 



direction (the left in FIG. |3J) when they are sheared by 
dislocations, the loops on the — y (the right) side of the 
precipitate cannot move across the precipitate to collide 
each other there. This phenomenon may be related to 
plastic cavitation which is often observed near precipi- 
tates 19], and also to the rotational plastic flow around 
a nondeformable particle poj . 

To conclude, the present study indicates that the dy- 
namics of dislocation-obstacle interaction is much richer 
than considered before. Cross-slip plays the crucial role 
in a bypass mechanism which is quite different to the 
Orowan mechanism. We have clarified the condition that 
determines which bypass mechanism occurs. Namely, 
the shear stress which is asymmetric with respect to the 
glide plane causes cross-slip of the screw dipole and the 
Hirsch mechanism follows. Further development along 
the line of the present study will be the relevance of the 
Hirsch mechanism to macroscopic plasticity. A three- 
dimensional continuum approach, which is referred to as 
"dislocation dynamics", is now developing to compute 
the many-body dynamics of dislocations and obstacles 
[22^ . It may be friutful that the present results are suit- 
ably incorporated to such continuum approaches. In ad- 
dition, it should be remarked that the dislocation dy- 
namics is closely related to the interface dynamics in the 
quenched disorder, which is extensively investigated in 
the field of nonlinear physics 23]. Dislocation physics 
would be much richer in communication with the differ- 
ent research field. 
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